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Abstract 

A four dimensional fermion determinant is presented as a path integral of 
the exponent of a local five dimensional action describing constrained bosonic 
system. The construction is carried out both in the continuum theory and in 
the lattice model. 


1 Introduction 


The problem of bosonization of fermionic theories was studied by many authors (see 
e.g. |]I|, PI, p[, PI, PI), mainly in the context of two-dimensional models. There 
were also some attempts to extend this procedure to higher dimensions [^, P, p|, 
l0| , but at present they lead only to partial success. 

The problem of bosonization is of a particular importance for lattice gauge theori¬ 
es as lattice QCD, because the fermionic degrees of freedom complicate numerical 
simulations enormously. 

Recently M.Lusher 0 proposed the algorithm for the approximate inversion of 
the QCD fermion determinant replacing it by an infinite series of bosonic determi¬ 
nants. In the present paper I will describe an alternative approach which allows 
to write the exact expression for the fermion determinant as a path integral of the 
exponent of a local bosonic action. In my approach a four dimensional fermionic 
system is replaced by a hve dimensional constrained bosonic one. 

In the second section I present the construction for the continuum theory. In this 
case the proof of equivalence will be given in the framework of perturbation theory 
as the problem of nonperturbative definition of a continuum path integral is at 
present beyond our possibilities. In the third section the corresponding procedure 
will be given for the lattice QCD where the equivalence may be proven at the 
nonperturbative level as well. 
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2 Perturbative bosonization in the continuum the¬ 


ory 


We consider the following fermion determinant 

det{D + m)'^= f exp{-^ f d‘^x[^l^{x){Dap + (1) 


2=1 ' 


where 




( 2 ) 


belongs to the Lie algebra of some compact gronp {SU{3) for QCD). Indices a 
denote spinorial and colonr components. The eq. ( describes n degenerate fermion 
flavours interacting vectorialy with the Yang-Mills held. To make contact with the 
discussion of lattice models we assume that the four dimensional space is Euclidean 
and the 7 ^ matrices {fi = 1,2...5) are Hermitean. As all the discussion in this 
section will be restricted to perturbation theory all the equations are understood 
as a series over g, and the Gaussian path integrals are dehned axiomatically as in 
ref. 1^, l]!^. Therefore we need not to bother about the convergence of the path 
integrals and all the conclusions are valid (perturbatively) for any number of fermion 
havours. However in the next section where we deal with the nonperturbative proof 
for lattice models we have to be more careful at this point and to make all the 
integrals convergent we restrict our discussion to the case of even number of havours. 
For simplicity we take n = 2. 

We present the determinant of the Dirac operator as the determinant of the 
Hermitean operator by using the identity 


Therefore 


det(D + m) = det[ 75 (D -|- m)] 


det{D + mY = det( 75 (D -|- = det(—-|- rn?) 


(3) 


(4) 


Being the square of the Hermitean operator the operator —D^ + rrA is positive 
dehnite. 

Let us introduce the following action depending on hve dimensional bosonic helds 

S = J dtd^x{X<i)*^{t,x)dt(j)a(t,x) + (j)*^(t,x)[-D^{x) + m\i 3 (j)/ 3 (t,x)} (5) 

Here A is a small positive parameter which eventually will be put equal to zero. The 
helds 0 Q, have the same spinorial and colour structure as the helds 'ifja- 

We claim that the determinant (|l|) (for n = 2) may be presented as the following 
path integral: 


det(D-|-m)^ = lirn y exp{—S'-|-i j dt j d^x(jY{x,t)Xa{x) 


+ 


nl 


J ^ xxl{x)(j)a{x,t)}d(j)ld(l)f^dxldxg 


( 6 ) 
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The fields x{^) ^ I'ol® of Lagrange multipliers imposing the constraints 


rl fl 

dt(j)*{t,x) = / dt(f){t,x) = 0 

D JO 


The Green functions of the fields 0 defined by the action(|[) are retarded: 

G(/c,p) = T ^ 


i\k + + m? 


G{t,p) = — / exp{ikt}G{k,p)dk] G{t) = 0,f < 0 
27r J 

Integrating in the eq.(^) over 0*, (f) one gets 

I = det(A(9t — + rn?)~^ x 


(7) 

( 8 ) 
(9) 

( 10 ) 


J / d‘^^d‘^yXaix)[^dt - D'^ + m‘^]J{x,y, s - t)xf3{y)}dx*dx 

Due to the fact that the Green functions of the fields (f> are retarded, for any A dif¬ 
ferent from zero the first factor in the eq.( IT) is equal to one. Indeed a perturbative 
expansion of this determinant generates the diagrams of the type 


n. 


dk 


{iXk + p\ + m?)... {iXk -|- -|- m 


( 11 ) 


where F is some polynomial over pj, k. Note that the fields A do not depend on s, 
therefore the diagrams generated by eq. m are the cycles with zero fifth component 
of momenta. The integral over k is convergent and all the poles are in the upper 
complex half plane. Glosing the integration contour in the lower half plane one sees 
that the integral is equal to zero. In fact one has to be slightly more careful at this 
point, as due to momentum conservation all these diagrams are multiplied by the 
singular constant factor ~ h(0). To give a precise meaning to this expression some 
infrared cut-off must be introduced to the eq. ( ^. In the next section we shall 
show that a finite lattice may serve this purpose. 

The second factor is not singular in the limit A —0 and one can calculate it by 
putting A = 0. Let us prove it in some more details. 

Expansion of the exponent in the eq. (|I0|) over g generates the terms of the form 



dk 


exp{ik{s - t)}P{Xk,pi ...pn) 
{iXk + Pi + rn?) ... {iXk + Pn + ’^^) 


( 12 ) 


where P is some polynomial over k,pi. Integrating over s and t one gets 


r (2 — e — e^^)P(Xk,pi) 

J {k — ieY{iXk + Pi+ m?)... {iXk + Pn + w?) 


The pole at A: = 0 is spurious. It is compensated by the zero of the nominator. 
So we are free to define the roundabout of this pole arbitrarily. We choose the 
prescription k k—ie. With this prescription the part of the integrand proportional 
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to (2 —exp{—z/c}) gives zero contribution as in this case we may close the integration 
contour in the lower half plane and all the poles are situated in the upper half plane. 
The last term, proportional to exp{i/i;} is calculated by closing the contour in the 
upper half plane. The integral is equal to the sum of the residues at the poles at 
= 0,/c = + m^). The contribution of the poles at /c = + m?) 

vanishes at A —0. The contribution of the pole at /c = 0 gives 


_ P{k,pi) _ 

{iXk + Pi + rn?) ... {iXk + Pn + m?) 

Therefore the integral m acquires a form 

1 = J exp{- J x*a{x)[-D^ + m\^{x,y)xi 3 {y)dxdy + O{X)}dx*adx 0 

Integrating over y*, y we get 


(14) 


(15) 


I = det(-D^ + m^) + 0(A) 


(16) 


Therefore 


lim / = det(Zi) + m)^ 

A^O 


(17) 


3 Bosonization on the lattice 

In this section we generalize the construction given above to the case of lattice 
models, where the proof of equivalence can be given without any references to per¬ 
turbation theory. 

We again present the determinant of the Dirac operator as the determinant of 
the Hermitean operator by using the identity 

det(.D-|-m) = det[ 75 (.D-|-m)], D = (18) 

In eq.( is the lattice covariant derivative 

= ^[U^{x)^p{x + a^) - Uf,{x)'tp{x - a^] (19) 

is a lattice gauge field. We consider a finite lattice with periodic boundary 
conditions. 

As it has been already mentioned, to provide the positivity of the determinant 
we consider the case of two degenerate fermion flavours interacting vectorialy with 
the Yang-Mills held. We shall not deal explicitely with the problem of fermion 
doubling. All the reasonings are trivially extended to the model improved by adding 
for example the Wilson term. 

It is convinient to present the fermion determinant in the following form 

det[75(£)-I-m)]^ = f exp{a‘^'^'ijj{x){D‘^ — m‘^)^jJ{x)}d^jJd^jJ (20) 

J rr 
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We shall prove that in analogy with the continuum case the integral over fermionic 
fields ^ can be replaced by an integral of a five dimensional lattice bosonic action. 
The spatial components x are defined as above. The fifth component t to be defined 
on the one dimensional lattice of the length L with the lattice spacing b: 

L = 2Nb, -N <n<N (21) 

We choose b in such a way that b « a and in the continuum limit 

2Nb‘^ = Lb^O ( 22 ) 


The key ingredient of the construction presented in the preceeding section was the 
using of retarded Green functions for the fields 0 providing triviality of the deter¬ 
minant in the eq.(W). This phenomenon is a continuum analog of the well known 
property of a triangular matrix: the determinant of a triangular matrix is equal to 
the product of the diagonal elements. Therefore to achieve the same goal in the 
lattice case we choose the lattice derivative with respect to the hfth coordinate t in 
the form of a triangular matrix. More precisely we are going to prove the following 
equality 

cpja^X !~ rri^)i^{x)}dipdilj = (23) 


exT 


X 


= lim 

A^0,b- 


exp{a'^6 


N 


n=-N+l X ^ 


0n(x) + 0*(x)(Z)^ - m^)0„(x)- 


-^iKix)xix) + x*ix)(pnix))]}d(j)*^d(j)ndx*dx. 

Here (j)n{x) and x{x) are bosonic helds which carry the same spinorial and colour 
indices as the helds 0(x). The helds (j)n{x) satisfy free boundary conditions i.e. 


0„ = 0, n < —N, n > N (24) 

The operator —D’^ + rn^ is Hermitean and can be diagonalized by a unitary trans¬ 
formation. As it does not depend on t this transformation will make the exponent 
in the r.h.s. of eq.( ^) diagonal with respect to all variables except for t. After this 
transformation the r.h.s. of eq.( acquires the form 

1= lim /(A, 6) 


/(A,4)= /exp{i 

„=_7V+1 a ® 

-^(0rx“+xr0“)]}#:d0ndx*dx; (t>N+i = o ( 25 ) 

Index a refers now to the eigenstates of the operator —D^ + being the 

corresponding eigenvalues. Obviously > 0. 
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The integral ( is convergent as the real part of the action in the exponent is 
positive. It is easy to see by rewriting the action in terms of Fourier components: 

1 — 

S = — [” {0*“(fc)[-A(e-'"' - l)r' + B'^]r{k)] + 

2 n 

0 

27rv L 

The real part is 

ReS = T {^*^{k)[-X{coskb - l)r^ + B‘^]^^{k)]dk (27) 

27T J-^ 

o 

As A < 0 and B°‘ > 0, ReS > 0. 

Performing in the eq.( |^) the integration over (l)n one gets 

r -b ^ 

I{X,b) = det{C exp{—J2 X*‘^{C'^)JnX'^}dx*dx ( 28 ) 

o n,m=-N+l 

Where Ca{k) is the kernel of the quadratic form in the eq.( 

As follows from the eq. ( [^) 

det((P) =ndet(C„) (29) 

a 

In the coordinate space Ca is a triangular matrix with the diagonal elements 

~ (A + Bab) (30) 

Therefore 

det(C') = exp{^ln(A + Bab)‘^^} (31) 

a 

Separating the constant term one gets 

det(C') = exp{2A^6A-^ E + 0{Lb)} = (32) 

= exp{X-^LTT [-D^ + m^] + 0{Lb)} 

Using the explicit form of the operator D one can easily verify that Tr [ —D^ + m^] 
is a nonessential constant. It follows also from the fact that the trace of a local 
operator is local. It has to be gauge invariant and a polynomial of the second order 
in the fields A^. The only possible solution is a constant. So we can include det(C'“^) 
into normalization constant. 

To get an explicit form of the remaining terms it is sufficient to find the stationary 
point of the exponent in the eq.( ^). The corresponding classical equations look as 
follows 

- A(C+i - <PT)b-^ + B^c^T + = 0; n^N (33) 

A(C - = 0; nj^-N + 1 
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- - iL-^x*"" = 0 (34) 

C7V+1^“^A + 0"jv+l^“ + = 0 

for small b eq.s( may be approximated by the differential equations 

- Xdt(j)*^ + + zL-^x*" = 0 (35) 

A<9t0“ + S>“ + zL-5x“ = 0 
whereas eq.s( |^) play the role of boundary conditions: 

) = 0; r(-\) = 0 (36) 


The solution of these eq.s is 


■(■“(t) = - <!xp{S“A-‘(t - ^)}) 


r(i) = 




^^__-X“(l-exp{-B“A-‘(i +j)}) 


(37) 


I ^ g„[^ “ exp{^(« - ^)}]dt}dx’dx = (38) 


Substituting these solutions to the eq. ( ^) we get in the limit 6 —0 

L 

lim/(A,6)= /exp{- p^ Y. 

= / exp{-i:x“(B“)"‘x”ll-;g^[l+<ixp{-B“A->L}]]}cix*<ix 

In the limit A —0 this equation reduces to 

/= lim /(A,6) =exp{-^x*“(i?“)-'x“W“^^x“ = det(-I)2 + m2) (39) 

A—>0,0—>0 


The equality (^^ is proven. The hnal recipe is the following: Any gauge invariant 
observable may be expressed as a purely bosonic path integral of the form 


Z = ^ \nn Jexp{a‘b f] + 

n=-V+l X 


\^0,b^0 

a2 2 


+K{^)iD - m )pn{.x) -^(0;(x)x(x) + x*{x)(j)n{x))] + 


Vl 


+L{U) + s.t.}d(j)'^d(j)ndx*dx 


(40) 


where L{U) is the lattice Yang-Mills action and s.t. stands for the source term 
depending on the helds All the helds are bosonic and the integration 

goes over the fields 0*, pn dehned on the one-dimensional lattice with length L and 
free boundary conditions. 
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4 Discussion 


In this paper we showed that a four dimensional fermion determinant can be written 
as a path integral of the exponent of a hve dimensional local bosonic action. In the 
same way one can present a two dimensional fermion determinant as a path integral 
for three-dimensional bosonic theory. In the case of lattice models this procedure 
leads to a well dehned bosonic path integral. No numerical simulations in this 
approach have been tried so far and it would be very important to see how the 
method works in practical calculations. In this case one need not of course to take 
the limit A —0, 6 —0, but the following conditions b « a, b << A, b << L~^ 
must be respected. 
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